Galilean-invariant scalar fields can strengthen gravitational lensing 
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The mystery of dark energy suggests that there is new gravitational physics on long length scales. 
Yet light degrees of freedom in gravity are strictly limited by Solar System observations. We can 
resolve this apparent contradiction by adding a Galilean-invariant scalar field to gravity. Called 
Galileons, these scalars have strong self-interactions near overdensities, like the Solar System, that 
suppress their dynamical effect. These nonlinearities are weak on cosmological scales, permitting 
new physics to operate. In this Letter, we point out that a massive gravity inspired coupling of 
Galileons to stress energy gravity can have a surprising consequence: enhanced gravitational lensing. 
Because the enhancement appears at a fixed scaled location for a wide range of dark matter halo 
masses, stacked cluster analysis of weak lensing data should be able to detect or constrain this effect. 



Our understanding of cosmology has been profoundly 
affected by the discovery of cosmological acceleration. It 
may signal a breakdown of General Relativity on long 
length scales. This has initiated a search for consistent 
modifications of GR. The leading models for modifying 
gravity are scalar-tensor theories: chameleonic / f(R) 
theories [U [2] and the Dvali-Gabadadze-Porrati (DGP) 
model [31 and its descendants [1]. Until now, models 
have assumed that the scalar couples only to the trace of 
the stress-energy tensor. Since radiation's stress-energy 
is trace-free, gravitational lensing is unaltered. 

The scalar field found in the decoupling limit of the 
DGP model, tt [5], has an intriguing quality: it is 
galilean invariant in the action. That is, the scalar 
part of the action is unchanged under the replacement 
TT — >■ TT + c + fopX^, where c and the &^ are arbitrary con- 
stants. This galilean symmetry can arise as a manifes- 
tation of higher-dimensional symmetries [5] , emerge as a 
consequence of giving the graviton a mass [THS] , or simply 
be posited as a foundation for model building [TU] . Fields 
with galilean-invariant actions are special: they are a 
symmetry-protected set of derivatively self-coupled fields 
with higher-order derivative actions, but with equations 
of motion that have only two derivatives operating on the 
field at a time. Equations of motions with more than two 
time derivatives are in danger of being ill-defined. Fields 
with this symmetry are broadly called galileons. Galileon 
equations of motion contain derivative terms raised to 
higher powers. The non-linearities introduced by these 
terms allow galileons to exhibit the Vainshtein mecha- 
nism [TT]: the scalar field becomes strongly coupled to 
itself near matter sources. This suppresses its gradients. 
Since the scalar force comes from gradients, their dynam- 
ical influence is suppressed near matter sources. This al- 
lows galileons to pass solar system tests while still having 
non-trivial effects on longer length scales. These effects 
include observables, like extra large scale structure and 
faster peculiar velocities [12]. 

In the decoupling limit of massive gravity [5] , de Rham 



et al. find a galileon-type theory with an additional cou- 
pling to stress-energy. This has a profound consequence: 
they are able to degravitate [3 [TB] , or suppress the back- 
ground curvature caused by, the cosmological constant 
at the linearized level. In this letter, we point out that 
couplings of the form described in [8] also have a strik- 
ing phenomenological consequence: they can significantly 
strengthen gravitational lensing relative to GR. 

The basic features of this enhancement are as follows. 
For a spherically symmetric source, it vanishes as r — >■ 0, 
giving negligible Parameterized Post Newtonian (PPN) 
effects. It also tends to zero as r — >■ oo, the limit where 
the dynamical effect of the field is largest. For the pa- 
rameters of the massive gravity model, the lensing shear 
is enhanced ~ 5% relative to GR for any spherically sym- 
metric mass configuration. The increased shear occurs at 
an intermediate length scale within the strong coupling 
radius of the theory, the so-called Vainshtein radius ~ see 
Fig. [l] This radius is given by = (?'s?'c)^^'^j where 
is the Schwarzschild radius of the source and Tc is the 
Compton wavelength associated with the graviton, typi- 
cally ~ c/Hq. For the sun, ^ kpc; for a typical galaxy 

~ Mpc; and for a galaxy cluster, ^ 10 Mpc. In 
the NFW profile, the change in shear is at the percent 
level for a wide range of radii (Fig. [2]). This lensing 
effect is qualitatively different from the parametrized de- 
viations from GR discussed in e.g. [E]: it is a locahzed, 
inherently nonlinear effect that disappears on long length 
scales and in linearized perturbation theory. It appears 
at length scales that are very well measured by galaxy 
surveys. The effect is nearly constant in r2oo units for dif- 
ferent halo masses and concentrations. Hence, it should 
be possible to discover or constrain this effect by stacked 
analysis of many halos' weak lensing data rescaled by 
their virial radii. Planned experiments like the Large 
Synoptic Survey Telescope (LSST) [H] should have suf- 
ficient depth to observe this effect. 

Enhancing the lensing potential: For the decou- 
pling limit galileon-type scalar field vr, called the helicity- 
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graviton, that arises in theories of massive gravity [HUHl) 
the couphng of the field to stress-energy has the form 



(1) 



where a and /3 are 0(1) dimensionless coefficients, A3 = 
{Mpim?gY/'^ is the strong couphng scale of the theory, 
Mpi = (l/G)i/2 is the Planck mass, and rug is the mass 
of the graviton. For our estimates we will take the gravi- 



ton to have a Compton wavelength, rc 



the Hubble scale today; and we will work in units where 
G ~ c — h — 1. This "Einstein frame" result is the 
simplest version of a class of theories studied in [S] . In 
this limit, the metric can be diagonalized and the scalar's 
effect more easily isolated. Despite these complications, 
it is clear from the derivations in |8l [9] that the metric 
whose geodesies determine the paths of photons is the 
one that includes both tensor (/i^i/) and scalar (9^7ri9^7r) 
parts. Earlier studies [Sllinj of galileon fields did not con- 
tain the coupling oc (Stt)^ although the absence of this 
coupling means that the tt field's stress-energy coupling 
is not obviously invariant under the galilean symmetry. 

As pointed out in [3], this novel coupling permits the 
degravitation of a small cosmological constant in the de- 
coupling limit. In this note, we point out that this cou- 
pling has another consequence: the enhancement of the 
gravitational lensing potential. 

For linearized GR, we have hoo = hij — ^Sij. For 
lensing in standard GR, the relevant potential is then 
given by $l = |(<I> — ^I'). In the presence of a spher- 
ically symmetric mass distribution, galileons generically 
have a non-trivial drTr and an approximately vanishing tt. 
The additional coupling changes the equations of motion 
slightly, but tt — >• is still a good solution. The extra 
coupling in the lagrangian implies that the potential $ 
is modified, leading to a fractional change TZ{r) in the 
lensing potential $l given by : 



$^$-t-A$, 7^(r) 



$l[GR] 



Ai 



(2) 

For our estimates, we will work with a general galileon 
theory [T0| , using coefficients consistent with [H [9] and 
including the extra stress-energy coupling found in [51 H] 
and given in Eqn. [T] The scalar part of this theory then 
has the lagrangian 



r 



(3) 



In this equation, we have abbreviated some expressions: 
((Stt)^ = d^ird^TT and n(j = d^d^n. We have also in- 
cluded five dimensionless 0(1) coefficients, a, /3, 77, /i 



Modification to point mass 3D lensing potential 
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FIG. 1: Radial dependence of the fractional change in the 
lensing potential, Eqn. [Sj for a point-like central mass, assum- 
ing ai = —1/2 and 02 = 1/2 (which gives a = 1, (5 = \, rj = \, 
fi = 3/2, and v = 1/2) in the scalar field equations. The ra- 
dius is scaled by r, = {vsr'^y^^, where Vs is the Schwarzschild 
radius of the source and rc is the Compton wavelength (or 
inverse mass) of the graviton, typically ~ c/Hq. For the sun, 
r, ~ kpc; for a typical galaxy r« ~ Mpc; and for a galaxy 
cluster, r» ~ 10 Mpc. The peak change of ~4% is achieved 
for r ~ 0.33 r. 



and v. Although gradients of n are suppressed near 
matter sources by the Vainshtein mechanism, the ap- 
pearance of the small scale A;^^ in the gradients' cou- 
pling to stress-energy permits A$ to become large. As 
we will see, the fractional change in $ is largest when 
tt' = dn/dr cx r^^/^. In spherical symmetry, the equa- 
tion of motion for tt becomes an algebraic equation for 
tt'. This equation is |10] 



3V 



+ 



aGM{r) 



(4) 

This admits a general closed form solution which is too 
lengthy to reproduce here. We have included 5 free co- 
efficients thus far, but in the massive gravity |8j case 
these are derived from just two parameters, oi and a2' 
a = — 2ai, P = 2a2, ?? = 4a^, = — 6aia2, and v = 2a\. 
(In [3], there is also a third free parameter, 03. When 
03 7^ 0, the action cannot be diagonalized into scalar and 
tensor components. Since this makes the physics more 
difficult to understand and is unnecessary to our pur- 
poses, we leave 03 = 0). This reduction of the parameter 
space gives a form of the solution for 7r'(r) that is differ- 
ent and simpler than the general cubic solution, due to 
a cancellation that occurs when 2/i^ = Q-qv. Note also 
that [To] finds general constraints on the parameters; for 
instance, ai < is required for radial perturbative sta- 
bility. We will specialize to the ai , 02 parameters for the 
remainder of this paper. The solution to Eqn. [4] as a 
fraction of the Newtonian force, 5*', is given in terms of 
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{2GMrlf/^, by 
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Next, we insert Eqn. [5] into Eqn. [2] and study its 
behavior for a point mass. The first thing to check is that 
the lensing modification vanishes near the origin, since 
gravitational lensing in this regime is tightly constrained 
by various PPN tests. We need 



dn 
dr 



(r - 0) oc r", 



n > 



(6) 



so that the behavior of the ratio A$/$l — >■ as r — )■ 0. 
This is what we find. Interestingly, the v = case - 
which recovers the galileon theory that emerges in the 
DGP model - has n — —1/2. This implies that the mod- 
ification to lensing from a DGP-like scalar would be non- 
zero at the origin. Since the enhancement amplitude is 
independent of A3, it persists even in the rUg — >■ 00 limit. 
This is forbidden by numerous PPN tests of GR. So in- 
clusion of higher-order terms in the galileon lagrangian 
was critical for finding an effect that is not already ruled 
out. This degree of non-linearity arises naturally in 
For our solution, the behavior near zero is given by 



dr 



(r ^ 0) = 



1/3 



r 
272 



a2 



), (7) 



i.e., approaching a small constant near r = 0, giving an 
n = scaling in Eqn. [6j Thus our solution does not 
violate solar system tests. 

The other limit to check is r — >■ 00. Here again, the 
ratio vanishes, since galileon theories generically recover 
7r'(r — 7^ 00) oc l/r^, so it scales as for large r. 

These limiting behaviors imply that the solution must 
at some point pass through the r~^^'^ scaling that will 
give a A$ with the same radial scaling as $l[GR] and 
hence a finite rescaling of the strength of gravitational 
lensing. For the parameters of the massive gravity model 
and general r, the fractional change in the 3D lensing 
potential is given by 
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where x — r/r.j,, = (2GMr^)^/'^. Note that the ratio 
takes a particularly simple form for ai — —1/2, a2 — 
1/2; we will make this choice in our plots. This choice 
also gives the same long-distance dynamics as the DGP 
model (i.e., 7r'(r)/*'^(r) -)■ 1/3 for r > r*). We have 
plotted TZ{r) in Fig. [T] 7?.(r) reaches a maximum at 
xo = rjr, = ((2\/3 - 3)a2/18a2)i/3 given by 



n{xo) 



1 

12 



(2a/3-3) ~ 0.04. 



(9) 



Modification to NFW tangential stiear 
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FIG. 2: The fractional change in the tangential shear for three 
NFW halo profiles as a function of radius scaled by the halo's 
virial radius (r/r2oo), assuming a\ — —1/2, 02 = 1/2 and 
rc — 3000 Mpc in the scalar field equations. These plots 
describe halos with any virial mass, M200 oc rfoQ. The three 
concentrations plotted are c = 8 (solid red), c = 6 (dashed 
blue), and c = 4 (dash-dotted black). 



This peak amplitude is independent of the parameter 
choices ai and 02 and can be regarded as a prediction of 
the theory. Though small, this modification gives a po- 
tentially observable modification to the tangential shear 
of extended halos; this is illustrated in Fig. [2] 

Weak lensing: The enhancement to lensing we are 
studying peaks on intermediate length scales. Weak lens- 
ing around galaxies and clusters is thus the best place 
to look for its effects. Hence, we calculate the effec- 
tive change in the tangential lensing shear caused by 
the galileon for a Navarro-Frenk- White (NFW) halo pro- 
file, following [TB]. We plot this for three different halo 
concentrations in Fig. [2] (N.B. Existing parameteriza- 
tions of modified gravity (e.g. jl5j) are designed to work 
on scales characterized by linear overdensities, fc < 0.1 
h/Mpc. The effects we are describing vanish on those 
scales, so they are not adequate to studying this effect.) 

For an NFW halo, the modification peaks at ^ 0.5 r2oo 
and at ~ 9r2oo for this parameter set, where rsoo is the 
virial radius. It depends quite weakly on halo concen- 
tration. Because M200 oc t^qq and oc M2QQ, the ef- 
fect peaks at the same locations, as measured in units 
of r2oo, for all M2oo- This makes the effect potentially 
observable: we can stack the lensing results from many 
clusters, scaled by their virial radii, and look for the effect 
to emerge statistically. The same reasoning also implies 
that the character of the modification will be redshift 
independent if the galileon's parameters do not depend 
strongly on cosmology. We should caution that this cos- 
mological behavior is not well understood. A simple es- 
timate of when the effect turns on is when the Universe 
comes within its own r*, which occurs around z ^ 1. So 
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our predictions are likely robust for 1 > z > 0. 

Detectability: Over the easily observable range r < 
1.5r2oo, (IA7I) - 1% (Fig. [2|. Taking this as a sig- 
nal above a known background, we can estimate what 
observations are needed to detect it. The GR shear 
at these radii is 7 ^ 10^^. Assuming a shape vari- 
ability of CT-y = 0.3, we find Nobs ~ 10'' observations 
are needed for S/N > 1. We can estimate Nobs — 
(A^ga//arcmin^)7Vjg^ggg^jgj^g. We can get Nobs ^ 10^ 
with an LSST-like depth of 40 galaxies / arcmin^ [T3] 
if we stack, e.g., > 5 x 10'* lenses that each subtend 5 
arcmin^. Our data cannot achieve this [T^- We are per- 
forming a more thorough study of detectability now |18j . 

Strong lensing: To see the effect of the galileon cou- 
pling on strong lensing, we can find a solution to Eqn. |4] 
for a singular isothermal sphere (SIS) and study its be- 
havior near r — 0, the strong lensing regime. The SIS has 
a density p{r) cx r^^ and a mass profile oc r. It turns out 
that the galileon-sourced 2D shear profile can be calcu- 
lated in closed form for the SIS in terms of hypergeomet- 
ric functions. For the SIS, the galileon-generated frac- 
tional increase in the lensing potential grows as {r/rc)^^^ 
for small r. This means that the galileon field generates 
an effective projected density profile S](^ ^ 0) cx C^"'^^'^, 
where ^ is the 2D radial distance from the center of the 
source after the linc-of-sight direction has been integrated 
out. Unfortunately, this component is quite small. For 
this additional source of effective surface density to gen- 
erate even a > 1% increase in the effective Einstein ra- 
dius, 9e, the mass per radius of the SIS would have to be 
> 10^"^AfQ/Mpc. This is unlikely to account for the ap- 
parent excess of lensing arcs seen in gravitational lensing 
surveys, e.g. [TO] . 

Conclusions: In this paper, we have given a first 
study of the modifications to gravitational lensing gen- 
erated by the inclusion of a new coupling of a scalar 
component of gravity to stress-energy. This coupling 
arises naturally in ghost-free theories of massive grav- 
ity [S], and is reasonable to include in phcnomcnological 
theories of galilean-invariant scalar fields. The generic 
effect of this coupling is to strengthen gravitational lens- 
ing on length scales ~ 0.5 r^,, where = (rgr'^Y^^; 
is the Schwarzschild radius of the source and is the 
Compton wavelength (or inverse mass) of the graviton, 
typically ~ c/Hq. For the sun, ~ kpc; for a typical 
galaxy ~ Mpc; and for a galaxy cluster, ~ 10 Mpc. 
The enhancement to tangential shear is at the percent 
level for the parameter combinations that appear in the 
massive graviton version of the galileon theory. The en- 
hancement appears at a fixed location in relation to a 
halo's virial radius for a wide range of masses and con- 
centrations. This should allow stacked analysis of weak 
lensing data to measure or constrain this effect. 
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